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ABSTHACT 


The effects of gravity- gradient torques during boom deployment 
maneuvers of a spinning spacecraft are examined. Two different 
configurations are considered: 

CD where the booms extend only along the hub principal axes. 

(2) where one or two booms are offset from the principal axes. 

For the special case of symmetric deployment (principal axes booms) 
the stability boundaries are determined and a stability chart is 
used to study the system behavior. Possible cases of instability 
during this type of maneuver are identified. In the second 
configuration an expression for gravity torque about the hub center 
of mass has been developed. The non-linear equations of motion 
are solved munerically and the substantial influence of the gravity 
torque during asymmetric deployment maneuvers is indicated. 
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NOMENCLATURE 


1*^2 


1 


(t) 
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3 
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V ' -X 


Time varying coefficients in the equations of motion 
for small out of plane Euler angles 

Offset of the control boom with end mass m^ from the 
d^, plane 

Unit vectors in the orbiting reference frame 

Unit vector along the geocentric position vector to 
the hub center of mass for asymmetric deployment 

Coefficient in the system characteristic equation 
for a rigid spinning spacecraft 

Offset of control boom with end mass m^^ from the 
d^.d^ plane 

Time varying coefficients in the equation for 
hfjhi for symmetric deployment 

Vector basis defined after the first Euler angle 
(Q^) rotation 

Coefficient in the system characteristic equation 
for a rigid spinning spacecraft 

fCl] 

IC 2 J = Constant vector appearing in the analytical 
solution for the Euler angles without gravity- 
gradient for symmetric deployment 

Vector basis defined after the second Euler angle 
(G^) rotation 

Extension rates along 1,2,3 principal axes, 
respectively 

Extension rate along 1,2 axes when inertia symmetry 
about spin axis is maintained and extension rate 
along all the three axes when they are equal 

Principal axes of the spacecraft 

Unit vectors along the body principal axes 



dm 


ElemencaL mass 


F = 

f(t) 

§0 

I 

K 

k 

L = 

\/Q 




M 


m 




® 1’“2 

N 


N^,N2,xM3 


B/2n^ 

C/JJ^ 

I(t)/I(t) 

Gravitational acceleration at the earths surface 

Angular momentum components in the body axes 

Instantanious values of principal moments of inertia 

Hub principal moments of inertia 

I^(t) = symmetric deployment 

(I^-I)/! = constant for rigid spinning spacecraft 

Gravitational constant for earth 

Angular momentum of hub about point Q 

An gular momentum of control mass with respect to 
point Q 

Boom lengths along the principal axes 1,2,3 
respectively 

Mass of main part of spacecraft 
Boom end mass 

Boom end masses along 1,2,3 principal axes 
respectively 

Control masses for asymmetric deployment 
Gravity torque 

Gravity-Gradient torque components 

^ 2 
2mc +2m_c- 



Constant appearing in the solution for angular 
momentum for torque free system 

Geocentric position vector to hub center of mass 

|r| » Radius of orbit for symmetric deployment 

Geocentric position vector of composite center of 
mass 

Position vector of control mass nu referenced to 
point Q 

Position vector of control mass m^ referenced to 
point Q- 

Position vector of composite center of mass 
referenced to point Q 

Position vector to elemental mass dm referenced to 
point Q 

Laplace Transform variable 
Time 

Coordinate of the control boom end mass m. along the 
dj^ axis 

Coordinate of the control boom and mass m. along the 
d^ axis 

(spin factor) 

Nutation angle 

Angular velocities about 1,2,3 axes resnectively 
(i = 1,2,3) 

Orbital angular rate 

Constant appearing in the solution for angular 
momentuBi for torque free system 
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D 


la 


( 0 ) 


“ Euler angles 
= 4mc^ 


(M+m^) /M+Sm) 


= ni2CM+iiLj^)/(M;+ZM) 

= -m^m^nrn-Zm) 


= Moment of inertia dyadic of satellite for symmetric 
deployraent 

=» Moment of inertia dyadic of hub for asymmetric 
deployment 

= Indicates differentiation with respect to t 
= Indicates initial conditions 


Subscripts 


Q 


reference point taken at hub center of mass 


I. INTRODUCTION 


A number of spacecraft have long telescoping appendages. These 
appendages might be on-board antennas which must be extended in 
orbit after the Initial injection sequence. The d3^amics of such 
spacecraft has been discussed in the recent literature in the 
absence of external disturbance torques.^ The purpose of this study 
is to determine the effects of gravity-gradient torques during the 
boom deplo3?ment maneuver. 

The first part of the current study will examine the effect of 

the gravity-gradient torque when the telescoping booms are deployed 

in pairs along the spacecraft principal axes. Possible use of such 

a deployment maneuver for detumbling a spacecraft has been examined 

2 

in a recent paper. From an application of Lyapunov's second 
method (using modified forms of the rotational kinetic energy as a 
Lyapunov function) sequences of boom extension maneuvers can be 

I 

determined so that the spacecraft will approach either of the two 
desired final states: close to a zero inertial angular velocity 
state, or a final spin rate about only one of the principal axes. 
This study did not consider the effects of external torques. 

For the special case of a gravitationally stabilized satellite 
librating in the orbit plane, the effect of gravity-gradient 

3 

during the deplo]nner.t maneuver has been studied previously. An 
approximate series solution has been obtained to simulate the dy- 
namics and the results compared with those of numerical integration. 



In the current study the non-linear equations of motion for a 
spinning spacecraft including the effects of gravity- grad lent 
torques are developed and these are solved numerically. For the 
special case of a symmetric spinning rigid spacecraft the stability 
chart previously developed can be used to study the system 
behavior.^ 

In the second configuration studied here, the system is 
assumed to consist of a central hub and one or two control masses 
offset from the principal axes. The dynamics of such spacecraft, 

6 — 

in torque free space, has been discussed in the recent literature. 
Reference 6 has examined the feasibility of a movable mass control 
device, for detumbling a large space station where a single internal 
mass is constrained to move along a linear track.. In Ref. 7 a 
control law for the boom mass position is obtained such that a 
quadratic cost functional involving the weighted components of the 
angular velocity plus the control is minimized when the final time 
is unspecified. In order to evaluate the gravity torque effects 
an expression for the torque based on asimilar procedure adopted in 
Ref. 8 is developed. The complete non-linear equations of motion 
with the gravity torque are obtained and the influence of the 
gravity torque due to asymmetry is illustrated. 


II, SYMMETRIC DEPLOYMENT 


A. Eulers Equations 

la the first configuration the booms extend along the principal 
axes of the spacecraft as shown in Fig. 1. It is assumed that the 
booms are massless and perfectly rigid. Using vector components in 
the d^,d^,d^ reference frame the Euler’s moment equations with 
time varying moments of inertia are 
&1 - (I2 - ^3) = ^3^ 

h2 - CI3 - “ N2 (1) 

S ■ “2“l ^^1 ■ ^2^ h 

where h^ = I^m^ (i - 1,2,3) • I^ = I^(t) are the principal moments 

of inertia, <»>. are the components of the inertial angular velocity 
i 

in the d^jd^.d^ frame and are the external torque components 

about the center of mass. In this analysis all the external 
torques of the system except the gravity-gradient torque are 
neglected so that N^ represent the gravity-gradient torque components 

B. Reference Frames 

Equations (1) will be expanded here in terms of the coordinates 
and unit vectors defined by Figs. 2. and 3., which together 
establish an orbiting reference frame and three attitude angles 

A A 

®1* ®2* ^3 body-fixed wiit vectors d^,d 2 »d 2 to the unit 

AAA 

vectors fixed in the orbiting reference frame a- .a^ja*. • Q, , ®3 



correspond to three successive positive rotations' about the vectors: 


A A 4^ 


A 


^1’ ^2’ *^3’ * Specifically, unit vector a-j^is directed 

along the radial line from the earth to the satellite mass center 0 
(local vertical) , unit vector a^ is directed along the trajectory 

A 

binormal (normal to the plane of orbit) and is defined to make 
a-, a„, a a right-handed orthogonal triad. We further assume 

X ^ J 

that the center of mass of the system moves in a circular orbit, so 

A 

that is along the path of the trajectory. 

The transformation from the principal body axes reference 
frame to the orbiting reference frame with the chosen 9^, Q^i ©2 
Euler angle sequence, becomes, after combining the three rotational 
matrices. 
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(2; 


where positive angles correspond to rotations in the positive 
right-hand sense as illustrated in Fig. 3. ’*s" represents the sine 

function and "c" represents the cosine function. 

By examination of Fig. 3, the expression for the inertial 
angular velocity, ilT, of the satellite in a circular orbit, can be 
written down as 

“ - V3 + 02b2 + + £^^3 (3 


where is the orbital angular rate. From consideration of the 


individual Euler angle rotations and Eq. (2) the following rela- 
tionships between the various unit vectors can be developed; 

A 

Cj - dj 

bj - ssii+cSjdj 

A A A A 

“ (S0^s02-c0^s02c02)d^+(s0^c0^+ce^s02s0^)d2+c0j^ce2d2 (■ 

After substitution of Eqs. (4) into Eq. (3) the inertial angular 
velocity components in the body axes can be expressed as follows: 

4 B * 

(jj^ = 02S02+£1 Cs0^s02-c0^s02C02)+9j^c02C02 
* • 

0^2 = 02C02+S2(s0j^c02+cOj^s02S02)-0^c02S02 

= 03-H2c03_c92-0^s02 ( 


C. Gravity-Gradient Torque 

The gravity-gradient torque N about the satellite mass center 
0 is given by^ 

^ 3k - n " ( 

N = — 3 a- X U • a K 

where 

with gg the gravitational acceleration at the earth’s surface and 

A 

Rq the radius of the earth, a^^ is the unit vector along the local 
vertical, D is the inertia dyadic about the satellite mass center 
and R is the orbit radius. For a satellite in a circular orbit. 






R 


* * '^If 2 s 

“3 ' “ ^ S ‘^3 ■'■ *^®3®®3^^’'3 


Equations (9) and (10) are the conq)lete non-linear equations of 
motion for the system. These are solved numerically as a system of 
six first order differential equations. 

E. Motion for Small Out of Plane Euler Angles 


For an initially symmetric satellite, if the mass symmetry 
(about the spin axis) is maintained during deployment, the 
equations of motion reduce to: 


h^+b(t)h2 = — b(t)I I^cG^se^sG^ 


h2-b(t)h^ = --^b(t)I I^cQ^sG^cO^ 


^3 “ ^0 ~ constant 


where 


I^(t) = l2Ct) = I(t) 


I--I 

t,(t) = _r — h 
^ ^ I-I "0 
3 


To study the stability of the syst*eia for small perturbations about 
the spin axis we examine the equations of motion assuming that 
Kl, l©ol are small (but not 9- which reflects the spin) . It 






should be noted that the angle between the normal to the orbit 


plane and the body 3 axis is a function of 0^^ and 0^ (see Fig. 3) . 
Under these assumptions, the expressions for the inertial angiilar 

velocity components, Eqs. (5), become: 

♦ 

= (02+fl0^)s02+(0j^-fi02)c02 

A. * 

“3 ^ ® 3 ^ 


and the corresponding components of inertial angular acceleration: 


• 


= (02-HI03^-030j_-K10293)s02+C03.-p02+Q302-^0j_03)c02 


«a ♦ 


0J2 = (02'H202“0303_'H20203) 003- (03^-^02+0302+^0^^03) sQ ^ 


“3 ” ®3 


(1 


Substitution of Eqs. (12), (13), (14) into Eq. (11) yields the 


eqtiations governing the imjtion in which |02 | remain small as 


cos ©3 + A 2 sin 0 ^ = 0 


A^^ sin 03 - A 2 cos ©3 = 0 


( 1 . 


where 


Aj_ = ©3_ - £n-03-b(t)}02 + f(t) (03-^02) + n03{03+b(t)} 

A, = 0’ + {n-0 -b(t)}0, + f(t)(0,+i2e,) + n0-{0-+b(t)} 



f(t) =■ lCt)/lCt:) 


In general these coupled equations cannot be solved in closed form. 

However the system behavior can be studied with the aid of the 

stability chart previously developed for a symmetrical spinning 
4 5 

satellite. ’ 

F. Stability Chart 

Since Eqs. (15) and (16) represent the equations of motion 
for a spinning spacecraft for which the principal moments of 
inertias are functions of time, these equations would also represent 
the equations of motion for a spinning spacecraft when the principal 
inertias are not functions of time, provided the terms due to time 
varying moments of inertia are modified so that: 

f(t) * 0 , b(t) = b =s ^ hjj = constant 

The modified equations are (under the requirement that Aj^=A 2 = 0 ) : 

0j^+(b+92“^)02+fi(b+02^®l = 0 

*’ * * * 

0--(b+e_-Q)G,+fl(b+0j0-+ iii- b0, = 0 (17) 

4 J 1 d 4 b}^ 4 

* 

Obseirve that 0^ is also a constant. We now define two important 
parameters, E, an inertia ratio and , a, a spin factor as follows; 
K= (I 3 -D/I 


a = 0_/£2 = (m -Q)/J 2 


(18) 



10 


With the aid of Eqs. (18), the parameter b can be written as 
b = Kto^ = KJ 2 (a+l) 

After introduction of these parameters Eqs. (17) become 
“ • 2 

0 ^+fl {K(a+l)+(a-l)} e^+fl {K(ct+l)+a} 0^ = 0 

0^-£2 {K(a+l)+(ci-l)} 0^-Hl^ {K(a+l)+a} 6^+3KS2^02 = 0 (19) 

If we take the Laplace Transforms of these two equations , the 
subsidiary equation can be arranged as follows t 


s^-H 2 ^ {K(a+l)+a} a{K(ori-l)+(a-l) }s 


ej^(s) 

-£ 2 {K(a+l)+(a-l)}s s^+fl^ {K(a+l)+a} + 




« 

s0^(O)+0^(O)‘Hl{K(a+l)+(o-l)}02(O) 


( 20 ) 


where Q^^Cs), © 2 ( 3 ) represent the Laplace Transforms of 0^^ and 02 , 
respectively, and ^^^(O), initial 

conditions . 

The stability of the system can be established by examining 
the characteristic equation associated with Eq. (20) which is of 
the form: 

- Bs^ + C = 0 


where 


= -[{K(a+l)+a} + f {K(cri-l) + (a-1)}" + 

^-r - {S(a+l)+a}^ + 3K{K(cri'l) + a} 

For stability, the roots of the characteristic equation 



Etust not have a positive real part. Two s 3 nnbols F^ and F^ are 
defined as; 



A brief study of Eq. (21) indicates that the condition under which 
there will be no positive real parts are; 

(1) F^ < 0 

(2) ?2 i" 

(3) (Fj^)^ - F^ > 0 

To obtain numerical restilts, the variation of the quantities 
2 

^1* ^2’ ^^1^ ” ^2 plotted in the a, K. plane. The boundary 

of the stable and unstable regions is established and the resists 
are illustrated in Fig. 4. The unstable regions are indicated by 
the cross-hatch marks. This stability chart has been previously 



The equations of motion for a spinning symmetric satellite 
with telescoping appendages in torque free space under the assump- 
tion that !o^l, I 02 I small are obtained from Eq,s,. (9) and (11) 
as: 


h2+b(t)h2 = 0 


(23. a) 

h2-b(t)h^ = 0 


(23.b) 

* ^1 
- ICt) 

^2 

®3 ■ lUT ®3 

(23, c] 

= i(t) 

^2 

®3 + UO ®3 

(23. d) 

0j - («3-a 


(23. e] 


The solution to the first two of the above equations has been 

1 


previously obtained as:' 


h^ = qQ* cos ( / b(t)dt + 


(24) 


h 2 = qQ* sin ( / b(t)dt + 


where qQ* and i|;q* are constants depending on the initial conditions 


For linear extension rates and under the assumption that the masses 
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emanate £i;om the center of mass of the spacecraft, the instantanious 
moments of inertia for symmetrical deployment become (see Fig. 1) : 


=• I(t) = I* + 2mc^t^ + ^m^c^^ t^ 
I^Ct) = + 4mc^t^ 

where 

\ = ^2 = ct, m^ - m^ = m 


(25) 


(26) 


^3 distances of the end masses m^^, m^, m^ from the 


satellite: mass center 0, respectively, c the extension rate along 
the 1, 2 symmetry axes, and c^ the extension rate along the 3 axis. 


Let 

p “ 4mc^ 

2 2 

p =» 2mc + 2m^C2 


(27) 


so that 


I(t) » I* + pt 


I^(t) = I3* + yt 


t t I (t)-I(t) 

f b(t)dt = f ^0 


( 28 ) 


Equation (24) can be written: 


h^Ct) = qQ* sin 



L*fvj 


r V 


where the constants qQ* and i|)q* are related to the initial condi- 
tions as follows: 


/hj' (0) + CO) 


-1 

' ‘>1^) 


After integration of Eq. ( 23 . c). 


©aCt) = tan‘^ [■^= 


flt + 03 (0) 


where 63(0) is the Initial condition at t = 0 . Equations ( 31 ) and 
( 29 ) are then substituted into Eqs. ( 23 . c) and ( 23 . d) with the 
result : 


0 ne = 

^ ^ i*+pt 


cosJ -~=- tan”^ h=V «t+ifi„*+0„(O) 


0 -K20 = — sin) tan“^ { ^ 

^ ^ i*+pt 


)- Qt+i|; *+0 (0) 


The variation of parameter technique for solving a system of 






first order differential equations is now applied to solve 
Eqs. (32). We seek a solution of the form 
^(t) * j[{^(t)u(t) 

where the ftmdamental matrix, jj(t) , is given by: 
jjr^(t) = cos flt sin J2t 

- sin Wt cos £2t 

and the vectors jiCt) and ^(t) are 
-1 

u(t) = / Jj (s) _^(s) ds + £ 

0 



In Eq. (35) , J.(t> represents the right-hand side of Eq. (32) and 
C is a constant vector. 
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4i = + G^CO) 


(37) 


The solution for 0^, 0^ can be obtained 'by first evaluating the 
integral in Eq. (36) and then substituting this result into 
Eq. (33). 



-siri(}> 


-COSi^ 


+ C, 


+ C, 


(38) 


^1* ^2 constants- which can be related to the initial conditions 
by 

” 0 ^( 0 ) 

(39) 

c, = ej(0) 

From consideration of Eq. (38), the following conclusions can 
be drawn; 

(1) For small t/Vl*/p , tan (t/l/l^) can be approximated by; 
ttfsFil , so that the frequency of oscillation is given by 
hQ/l* = I^*(jj2C0)/I*. This indicates that for high initial spin 
rates oscillations will be characterized by a high frequency mode. 


(2) For t sufficiently large, tan - can be approximated 

VPp Vi*/p 

^0 IT 

by a constant: -r . 3^ and 0- are then periodic with a fre- 

Vi^ ^ 

quency of J2 (the orbital frequency) . 

The analytical solution developed is plotted in Fig. 5. a and 
is compared with numerical integration including the effect of 
gravity-gradient, for an initial spin rate of 3 rad/sec. The 
behavior obseirved is a confirmation of the above two conclusions. 

As time increases the analytical solution becomes less valid due 
to the accumulated effects of the gravity-gradient torques. 

H. Numerical Results 

The non-linear equations of motion (9), (10) are solved 

numerically. The numerical integration is carried out using a 

9 

NOVA-840 computer, with RKGS subroutine. The subroutine RKGS 
solves the initial-value problem fay means of a fourth-order 
Runge-Kutta formula using the modification due to Gill. The 
integration procedure is stable and self-starting; that is, only 
the functional value at a single previous point is required to 
obtain the functional values ahead. For this reason it is easy to 
change the step site at any step in the calculation. The entire 
input of the procedure is: (1) lower and upper bound of the 

integration interval, initial increment of the independent 
variable, upper bound for the local tnincation error; (2) initial 
values for the dependent variable and weights for local truncation 


errors in each, component of the dependent variables; (3) the number 
of differential eqr itions in the system; (4) as external subroutine 
sub-programs, the computation of the right-hand side of the system 
of differential equations; for flexibility in output an output 
subroutine. 

To study the system for which » l®2l small the 

stability chart developed (Fig. 4) can be employed to analyze the 
following cases of interest; 

(1) Extension only along the '3' principal axis ; A previous 

2 

study in the application of telescoping booms for detumbling a 
spacecraft, to achieve a desired spin about one of the principal 
axes has suggested (based on a modified form of kinetic energy as 
a Lyapunov function) extension of booms along all three principal 
axes until the final desired spin rate is reached and then con- 
tinuing the extension of the set of booms along the nominal spin 

axis until the transverse components of the angular velocity , 

2 

reach an acceptably small amplitude. For extension of booms along 
the spin axis only and for a symmetrical spacecraft it is clear 
that in the stability chart the spin parameter, a, will remain 
constant during deployment, while the moment of inertia ratio, K, 
decreases. The value of K at any instant of time for this case is 
given by 


(I t^) 



where starred quantities represent the initial values and 



m = end mass 
c = extension rate 

From consideration of Eq, (40) as t becomes large, K, tends towards 

-1. In the upper-half of the stability chart, the stability 

boundary, F 2 = 0, tends toward K = -1*0 asymptotically. This 

implies that for a wide range of spin rates commonly used, extension 

only along the *3' axis could result in the system eventually 

being driven to an unstable state. Typical examples of such 

deployment maneuvers have been simulated and the responses shown in 

Figs. 5.b,c; 6.a,b; 7.a,b; 9.a,b; and 10. a, b. 

In Fig. 5.b the response of the and components of the 

angular velocity with and without the gravity-gradient torque are 

compared. The initial conditions, end masses, extension rate and 

» ! 

initial principal inertias are indicated in the figure. For' this 

case an initial (high) spin rata of 3,0 rad/sec is assumed. 

Clearly, the envelope of and with gravity-gradient torque 

present shows the growth in amplitude after 250 secs. For this 

case the response of the nutation angle (Fig, 5,c) indicates the 

potential instability during extension after 250 seconds. Since 
3 

a * 2.9 X 10 here, this case can not be indicated in the regions 
of the stability chart as plotted. However use of Eq. (40) shows 
that for t > 12 secs, the value of K becomes less than -0.8. 



Figures 6.a,b and 7.a,b illustrate the responses of 
and nutation angle for two different extension rates; 
ft/sec and =* 0i5 ft/sec, respectively, both with the same 
initial conditions. In either case, a = 1.0, and is a constant 
during the extension maneuver. Initially K =® 1.0. Figs. 6. a and 
7. a show that with the gravitational torque included, the com- 
ponent increases in amplitude after an initial tendency to decrease. 
After 60 secs for- the first case (Fig. 6.b) and 250 secs for the 
second case (Fig. 7.b), the nutation angle increases and almost 
doubles its initial value within 280 secs and 500 sacs, 
respectively. The time history of the slower extension case 
(0.5 ft/sec) is indicated on a redrawn version of the stability 
chart in Fig. 8. In this chart the horizontal dotted line parallel 
to the K axis passing through o * 1.0 corresponds to the 0.5 ft/sec 
extension rate discussed earlier in Figs. 7,a,b. The system 
crossed the stability boundary at t = 39 secs, which corresponds to 
approximately 20 feet of boom extension. It can be concluded from 
Fig. 7,b that’ although the system crosses the stability boundary 
at t =» 39 secs, the instability is not exhibited in the form of 
the growth of the nutation angle during the extension maneuver, 
until about 250 secs , when the nutation angle begins to grow in 
an exponential fashion. 

Two more cases of extension only along the *3' axis were 
simulated, beginning in the unstable region in the right-half of 


the stability chart, that is the unstable region for positive K. 

The responses for the two cases with extension rates of 4 ft/sec 
and 0.5 ft/sec are shown in Figs. 9.a,b and Figs. 10. a, b respec- 
tively. The responses are similar to the cases of Figs. 6 and 7. 
Even for the slower extension rate of 0.5 ft/sec, no appreciable 
increase in the nutation angle was observed as the system passes 
through the unstable region for K > 0. 

(2) Extension maneuver along all three axes 
2 

A previous study on detumbling a randomly spinning spacecraft 
using telescoping appendages had considered zero inertial angular 
velocity as a final desired state of the system. The authors 
using the rotational kinetic energy as a Lyapunov function had 
concluded that the necessary conditions for asymptotic stability in 
torque free space are satisfied for positive constant boom extension 
rates. When the dynamics for symmetric extension with the 
gravity-gradient torque is considered, it is observed that the 
parameter a ° ((a^-Q)/Q eventually tends towards -1,0 for suf- 
ficiently large time. If the same end masses and extension rates 

were employed along all three principal axes, the inertia ratio K 

✓ 

at any instant of time would be given by 
I *-I* 

K = -J (41) 

I*+pt 

2 . 

ji = 4m c ; m = end mass; c = extension rate 


where 


Equation (41) suggests that for sufficiently large time, K tends 
towards 0. Although the point (0, -1) in the stability chart is a 
bovmdary point, in practice, extension along all three axes with 
the same end masses and extension rates could never result in the 
system actually reaching this point. 

Typical examples of such deployment maneuvers have been 
simulated and the responses shown in Figs. 11 and 12. a, b. In 
Fig. 11 the system is initially in the stable region with 
o(0) =5.0 and K(0) = 0,6. The initial conditions are indicated 
in the figure. For this case the responses of the components of 
the angular velocity and the nutation angle show negligible dif- 
ference with gravity-gradient present or absent. The nutation 
angle is observed to be a constant. The time history for this case 
has been indicated in Fig. 8. by the dash-double dot curve 
beginning at the point (0.6,5). With the assumed extension rate of 
4 ft/sec and end mass of 0.01 slugs, the system moves rapidly down 
the stability chart. It enters the unstable region for positive K, 
at about 3 secs which corresponds to 32. feet of boom length. 

However it stays only very briefly in this unstable region 
reentering the stable region again in about 10 seconds. Clearly 
any such deployment strategy should ensure that the extension of 
booms is not terminated in this unstable region. 

Figs. 12. a, b show the response of a system initially in the 
unstable region with a(0) = 5.0 and K(0) = -0.6. Here a small 



difference in the response of was observed with gravity- 

gradient present and then absent. Within 630 secs of response time 
this difference is not noticeable within the scale shown in the 
figure. However the response of the nutation angle (Fig. 12. b) 
shows a slight growth for the case with gravity-gradient present 
after about 430 seconds. The time history for this case is 
indicated by the dash— dot cuirve beginning at the point (-0.6,5) in 
Fig. 8. It is observed that the system tends to the point (0,-1) 
from the left. The system briefly passes through the narrow strip 
of the stable region close to the K axis in the left-half of the 
stability chart, although never entering the stable region elsewhere 
even for the slow extension rate of 0.1 ft/sec used. 

(3) Extension along the 1,2 principal axes 

A number of spin-stabilized satellites have long appendages in 
the plane of rotation. Hughes^^ has studied the dynamics of the 
satellite during the deployment maneuver assuming torque free condi- 
tions. For symmetric extension along the 1,2 principal axis only 
the spin parameter, a = tends towards -1.0 for 

sufficiently large time. The inertia factor K at any instant of 
time during the deployment maneuver is given by 


I,*+yt^-I*- ^ t^ 
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ji = 4 me ; m = end mass; c = extension rate. 

From Eq. (42) it can be concluded that, given sufficiently large 
time the factor K tends towards 1.0. So one can expect that under 
suitable conditions extension only along the 1,2 principal axes 
would eventually drive the system to the point (1,-1), which lies 
in the unstable region for positive K in the stability chart 
(see Fig. 4). 

Two typical extension maneuvers for this case have been 
simulated and the responses presented in Figs. 13. a, b and 14, a, b. 
Fig 13, a shows the response of the transverse components of the 
angular velocity, for a system initially in the stable region; 

K(0) = 0,4 and a(0) =» 5.0. The initial conditions are indicated 
in the figure. It is seen that the component for the case with 
gravity-gradient torque, begins to increase after about 270 secs. 
The nutation angle (Fig. 13. b) likewise begins an exponential type 
growth at about the same time. The spin, is not affected since 
the gravity-gradient torque about the spin axis vanishes for 
symmetric deployment and the responses with or without the gravity- 
gradient torque are identical. The time history trajectory for 
this case has been indicated in Fig. 15 by the dash-double cross 
curve beginning at the point (0,4,5). For the assumed extension 
rate of 0.5 ft/sec and an end mass of 0.01 slugs the system crosses 
into the unstable region for positive K after 87 secs which 
corresponds to about 44 feet of extended boom. Although the 


system enters the unstable region after 87 secs, instability is not 
exhibited in the form of growth in the nutation angle until about 
270 seconds . 


Another case of deployment along the 1,2 axes beginning in 
the unstable region (a(0) = 5.0, K(0) = -0.6) has been simulated 
and the results presented in Figs. 14. a, b. A slow extension rate 
of 0.1 ft/sec has been assumed for this case. Fig. 14. a shows 
that with the gravity-gradient present the component becomes 
positive at about 500 secs and continues to grow further, whereas 
the response without gravity-gradient torque shows that tends 
towards zero. Fig. 14. b shows that with the gravity-gradient 
torque the nutation angle begins to grow after 360 seconds. The 
time history for this case is indicated in Fig. 15. by the dash-cross 
curve beginning at the point (-0.6,5). The system enters the 
stable region at t = 74 secs, leaves it again and reenters the 
unstable region for positive K at t = 179 secs. Although the 
system is moving through the unstable region for t > 179 secs, it 
takes about 360 secs before the instability is indicated as a 
pronounced growth in the nutation angle. 



III. ASYMMETRIC DEPLOYMENT 


A. Configuration, 


The system is assumed to consist of a rigid central htib 

(Fig. 16. a) with center of mass at point Q and one or two extendible 

telescoping booms with end masses m^ and m 2 , respectively- The mass 

along the boom lengths is assumed negligible in comparison with the 

end masses. It is assumed that, in general, the two booms will be 

offset from the hub principal axes with the coordinates a, b, c, d 

indicating the amount of offset. Previous studies have considered 

6 7 

this type of configuration for detiaabling a spacecraft. * It has 
been pointed out that for three axis optimal control more than one 
offset boom (orthogonal to each other) is required and that for two 
axis optimal control a single offset boom is sufficient.^ 

B. Development of Gravity Torqije Components. 


The general configuration of the two mass offset system is 
shown in Fig. 16. b. Whenever there is an asymmetric (internal) 
mass motion in a spacecraft system, the position of the composite 
center of mass and the orientation of the system principal axes 
will vary with time. The choice of the composite center of mass 
of the system as the reference point in the body, leads to time 
varying moments of inertia in the rotational equations of motion. 
This problem can be circumvented by choosing the vehicle (hub) 
center of mass as the reference point and the hub principal axes 
as the body reference frame. In Fig. 16. b, dj^,d 2 ,d^ represent the 


hub principal axes; Q, the hub canter of mass; R, the geocentric 

position vector of the hub center of mass; R the geocentric posi- 

tion vector of the composite center of mass, C; and ^ 2 ^ are the 

control masses whose position vectors relative to point Q are 

indicated by r. and r„, respectively, r is the position vector 

of the composite center of mass referenced to point Q. 

If r represents the vector from point Q to an elemental mass 

dm, then the torque applied by gravity about the hub center of 
. 8 

mass is 



k(R 4- r) 

— — 3 

R + r 


dm 


(43) 


2 

where k = SqRq with g^ the gravitational acceleration at the 
earth’s surface and 8^ the radius of the earth. Eq. (43) can be 
expressed as: 



(44) 


where R = |Rj and the Integration is taken over the total mass. 

Expanding the denominator of the integrand in a binomial series 

, , 8 
provides the approximation 


in which terms above the first degree in (R • r/R^) have been 
neglected in comparison with first degree terms. This approximation 
is justified by the observation that the satellite dimensions are 
small in comparison with the orbital radius. The main body can 
be treated as a continuum and the offset masses created as discrete 
particles so that the integral in Eq. (45) can be divided into an 
integral over the main body and a summation for Che end masses. 
Equation (45) can then be written as: 



body 


1 



dm + ~ 


3 


R 


2 

X 2 

i«l 


r. 

i 


1 


r2 


(46) 


By definition of the center of mass of the hub /rdm is zero. 
Eq. (46) then becomes 


"q j^3 1 


m 


* k * ^ — 
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r'i 


(47) 


where 

A. 

a£ - unit vector along the geocentric position vector of the 
main body center of mass, 

R - distance of the main body center of mass from the center 


of the earth. 



Li - moment of inertia dyadic of the main body, 
m 

Here it is important to observe that in general the factor Z m.r. 

i ^ 

does not vanish. Eqtiation (47) is an expression for tire gravity 
torque about the hub center of mass for a system with two offset 
control masses- Here it will be expanded to yield the gravity 
torque components in the hub reference frame for the specific 
configuration shown in Fig. 16. a. 

A 

We now make the assumptions that the unit vector a£ defined 
above can be considered equivalent to the unit vector along the 
composite system local vertical and that the distance, R, as 
defined above is equivalent to the orbital radius. These 
asstmiptlons are justified by the fact that for the choice of 
control masses as one percent of the satellite mass,^ the dis- 
placement between the htib center, 'of mass and the composite center 
of mass would be extremely small in comparison with the orbit 
radius. It is then possible to make use of the same set of Euler 
angles (Fig. 3) as defined for the symmetric deployment in 
Section II, now relating the orbiting reference frame (Fig. 2) to 
the hub principal axes reference frame, d^^jd^jd^. The trans- 
formation between these two reference frames is given by Eq. (2) 

A 

and the vector a£ is obtained from Eq. (4) and is given as follows 

A A A A 

aT * c0-c0-d,-c0oS0-do+s0_d, 


(4 



For the two mass system shown In Fig. 15 the control mass position 

vectors are given by; 

^ ^ ^ 

r^ =* a + b + z d^ 

_ * -s .. ( 49 : 

r^ =» X dj^ + c d^ + d d3 

Equations ( 48 ) and ( 49 ) are svibstituted into Eq. ( 47 ), resulting in 
the following gravity torque components 

Ni = ^ (I2-I3) C02S02S03 + ^ [m^ {-zc02S03-bse2> 

E E 

3 k 

+m2 {~dc02S03-cs02}] 2 

E 

+2382) (-zc02S03”bs02) + m2(xc02C03-cc02S03 

+ds 02 ) (-dc 02 S 03 -cs 02 ) ] ( 50 . a 

51c Ic 

N2 = -^ (1^-13) 002802^03 + ^ [m^^ {-2c02C03+as02} 

R R 

3 k 

-hii2 {-dc 02 C 034 xs 02 }] 3 [m^(ac02C03-bc02S03 

E 

+2382) (as02-zc02C03) + m2 (xc02c0 3-0062303 


4 ds 0 ») (xs0_-dc0„c0_) ] 


( 50 . b 



® (I 1 -I 2 ) (c®2^ + — 2 t®^(t)c92C02+ac02S02} 

R R 

3k 

-Hfl2 {cc62a02'^‘^®2®®3^^ 3 

R 

+zsQ^) (bc0^cd^+ac9^sQ^) +m^(xc6^c6^-cc9^sQ^ 

+ds02) (cc© 2C02'^‘^®2®®3^ ^ (50. c) 


C. Equations o£ Motion 

The complete equations of motion with telescoping type control 

booms in the presence of gravity torque are developed. The torque 

free equations of motion for the system have been previously 
7 

developed and this development is briefly reproduced here, after 
modification to include the gravity torque. 

The generalized vector equation of motion for such a system 
containing a central hub and moving connected masses can be 
written 

’ 1 

. a 


' '■q ”1 ^ 




(51) 


where refers to the gravity torque, Q refers to the reference 
point which is assumed to be at the center of mass of the hub, R 
is the inertial acceleration of the reference point and Is the 

position vector of mass, m^, with respect to point Q (Fig. 16. b.) 

It should also be noted that is the position vector of the com- 
posite system center of mass whose position will change with 
the movement of m^^ and m 2 > The composite c.m. is assumed to move 


In a circular orbit, and it is assumed that coxipling between orbital 
(translational) motion and the attitude dynamics is a higher order 
effect. 

Ihe angular momentum of the system measured with respect to 
point, Q, has three con 5 )onents , 


S/Q ^hLj_/Q ^m^/Q 


(52) 


where describes the momenttm of the hub, and describes 


the momentum of mass m^. The hub momentum may be expressed in 


terms of the hiib principal moments of inertia and angular velocity 
components as: 


H/q “ Wi 


(33) 


A A 


where d^, d^, d^ are unit vectors along the hub principal axes, 
and 


^m /Q “ “i ^""i i i C5A) 

where r^ describes the position of m^ relative to 

We will now consider the inertial acceleration of the 
reference point (?ig. 16) . 

t = R - r (55) 

c c 

I&ider the assumption that coupling between translational and 
rotational motion can be neglected, R^ = 0, and 
R = - r 

c 

From the definition of the system center of mass we can relate 


( 56 ) 



R =3 ** 

M + Em 

where M represents the hub mass and 
Em = 

After substituting Eqs. (53), (54) and (57) into Eq. (51) the 
following rotational equation results;^ 

^/q * “l == “2 ''2 * '2> 

where 

^1 * ’’^l 

Eq. (58) is then expanded using the familiar relationship, 

i _ _ 

^/Q dt body ^ H»/Q 

and for the specific geometry of Fig. 16. a, 

rfV A A 

r^ = a dj^ + b d^ + z d^ 
t2 = X di + c d^ 4- d ds 

The acceleration terms r. (i = 1,2) may be calculated by using 



= w X (u X r^) + u) X + 2uj X 

^\hody 

together with Eqs. ( 60 ) and ( 61 ) . 

The complete nonlinear equations of motion are obtained by 
expansion of Eq. ( 58 ) and substitution for the gravity torque 
components from Eqs. ( 50 . a), ( 50 . b), and ( 50 . c). 

• 2 2 ’ 

+ (I3-I2) <fl2“3 Kb +z ) - abm^ - azm^ - azw^di 

2 2 * 2 2 " 

+(b -z ) (*>2“3 *** *** ^‘*'3 ” “2 ^ bz] 

+U2 [c 4d ) -csoo^ - dxM^ - dxm^m^ + (c -d ) W2**^3 

. • 2 2 

+0x^2^013 - 2cxLt>2 - 2dxw3 ^“3 ”‘^2 "** V^3U(bc4dz) 

* • 

-(ac+bx) - (ad+xz) CU3 - (ad+zx) ^^^(03 + 2(bc-dz) 01^013 

* * 2 2 

+(ac+bx) 2dztii^ - 2bxn)2 - 2xziii3 + (bd+cz) (m^ ) 

+ cz] = ^ ^2-13)^9230239^ t®l Kzc02S03-bs02> 

R R 

3k 

+m2 {-dc92s92-cs02}] 3 [m^ (ac02cO3-bc02S03 

R 

+3302) (-zc02s93-bs02) + ^2 (XC02C03-CCO2S03 


+ds9„) (-dc0„s9«-cs0„) ] 


1^102 + " (a +2 ) W2 + bztii^-bzw^iu 

2 2 • 2 2*’ 

+abtii2W2 + Ca -2 ) tu^w^ - 2zzoi2 + azCu^ ) + az] - \i2^[ 

-(d^+x^)w2 + cdtu^ - cdm^(Ai2 + cxiii2t02 ■*" 

* 2 2** ^ ^ • 
-2x32112 + dxCoj^ “ 2(a3cfdz)Qi2 

+ (bd+cz)ci)2 ~ (bd'K:z)<i3j^tii2 + (acrHix) 012(1)2 2(ax-dz)o)2W2 

^ * 2 2*’** 

-2(a^dz)oi2 + (ad+xz) (oi^ + xz-zx] 

= ^ ■*■ \ [®2{-zc02c 93 + as02> 

3k 

+ni2 {“dc02C03 + xs02)] 2 

R 

+zs02> (as02-zc02C03) + m2(xc02C03-cc02S03 

' / 

+ds02) (xs02-dc02C03) ] 

I3W3 + <l2-I]_)(i)2tU2 “ U^Cazw^ + bzoJ2 - (a^+b^)o)3 

2 2 • " / 2 
+(b -a ” ^2(1)20)3 + bzto^Wjj^ + 2azoi^ + 2 bzbi 2 ^ ab((02 * 

-P2[<i3au2 + “ Cc^+x^)o)3 + (c^-x^)tU2«D2 “ <i3aj2“3 

+cdoi3W^ - 23001)3 + cxCw^ - 0)2 ) + ex] - P3[ (ad+xz )(iij^ 

+(bd+cz)oj2 - 2(bc+ax)ii)3+2 (bc-ax)oi2W2 " Cad+xz)u)2U3 

• • • 2 

+ (bd+cz)aj3ai^ + Iszvi^ + 2czoi2 “ 2331013 + (ac+bx) (oi^ ~<^2 

+ bx] = -^ (c 02 )^c© 3S03 + “ tnLj_{bc02c93 

R R 

3k 

+ ac02S03) + m2 {0002^03 + XC92S93}] 3 [m^ (3003003 

R 

- bc02S03 + 2303) (bo02C03 + 3062303) + m2(xc02O03 

- CC02S03 + ds02) (0003003 + 3CC02S03) ] 
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Equations (9) , which relate the Euler angular rates to the angular 
velocity components, together with the above Eqs. (63, a), (63. b) 
and (63. c) are the complete non-linear equations of motion for the 
two mass asymmetric system. Due to their complexity, no attempt 
has been made to obtain an analytical solution. 

D. Numerical Results 

The equations of motion developed are solved numerically 
using the EKGS subroutine outlined earlier in section II. In 
Ref. 7 a control law has been obtained for the boom end mass 
position such that a quadratic cost functional involving the 
weighted components of angular velocity plus the control is 
min imi zed when the final time is unspecified, assuming torque free 
conditions . 

*■ 

As an illustrative example the system parameters and Initial 
conditions are selected from Ref. 7 for a large space station. It 
is to be noted that the satellite mass and the control mass here 
are 32 times larger when conq>ared with, those used .in Ref. 6. 

= I 2 = I = 10.5 X 10^ slug-ft^ (1.42 X 10\g-m2) 

I^ * 15 X 10^ slug-ft^ (2.03 X lo\g-m^) 

M = 1.37 X 10^ slugs (6.21 x 10^ kg) 

= 1800 slugs (26112 Ibm) 

a = 65 ft (19-8m), b = 0 ft (Om) 

0)^(0) = 0.391 rad/sec, 0 ^ 2 ( 0 ) = 0.0 rad/sec, m^(0) = 0.314 rad/sec. 
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In all the cases in this study only a single boom (Z boom) offset 
from the spin axis has been considered. A specific example has 
been chosen from Ref. 7 for which the Z boom control law is stated 
as follows: 

z + (0.314)^ z = - 0.314 X 17.72 {kc^u^+kc^m^ } (64) 


where kc. and kc are constants chosen based on optimal 
7 

control theory. In this case: 


kCj^ = - 3.92, kc^ = 9.56. 


The solution to Eq. (64) with the initial conditions z(0) = 0, 

A 

z(0) = 0 for the torque free case is given in Ref. 7 as: 
z(t> = 17.72 {0.655 sin 0.314t + 0.417 cos 0.314t 

(65) 

_^-0. 54x0. 314t (^^ 417 0.428x0. 314 t+2. 056 sin 0. 428x0. 314t)} 

In the first part of the nvimerical study the effectiveness of 
the above control law in the presence of gravity torque is deter- 
mined. Two cases of interest have been identified and the system 
behavior for these two cases are compared with the behavior for 
the torque free system. 

Case.l; where the control law as given by Eq. (64) employs the 
actual components of the angular velocity tiii, 
presence of gravity torque. 

Case. 2: where the Z boom motion is according to the idealized 

control, the controller using the angular velocity components for 


the torque free system. For this hypothetical case the Z boom 
motion is input into the program as given by Eq. (65). 

Figure 17. a shows the Z boom motion for Case.l by the dash- 
cross curve, and for Case. 2 by the solid line. It is seen from 
Eq. (65) that the boom will experience a steady state oscillation 
after the initial transient for Case. 2. Observe that for this 
case within 10-15 seconds the transient part of this motion is 
removed leaving a remaining steady state oscillation. For Case.l 
the Z boom motion amplitude increases after 40 secs and would 
reach displacements of nearly 300 feet in 60 seconds (assuming 
that much boom length could be extended) . 

Fig. 17 .b shows the response of the component of the 
angular velocity. When gravity torque is neglected the amplit\ide 
of 0 )^ decreases with time and becomes close to zero in about 20 
seconds of response time, whereas for Case.l. and 2. the response 
shows a markedly different variation with becoming as high as 
-0.12 rad/sec and -1.5 rad/sec respectively within 60 seconds of 
response. Figs. 17. c shows a similar behavior for the 1^2 com- 
ponent. For the case without the gravity torque Im^l becomes 
close to zero as opposed to the large amplitude motion for Case.l 
and Case. 2. Figure 17. d indicates that is constant when 
gravity torque is neglected, whereas for Case. 2 the motion is 
periodic over a mean value of 0.6 rad/sec. For Case.l, correspond- 
ing to the actual Z boom motion, the response is somewhat 
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similar to Case. 2 until about 40 seconds after which it sxiddenly 
increases rapidly. This sudden Increase corresponds to the high 
amplitudes of the Z boom motion. Fig, 17. e. shows the response of 
the nutation angle. For all cases there is an initial tendency to 
reduce this angle. However after 8 secs, with gravity-gradient present 
there is a definite tendency to exceed the Initial value. 

For a single boom offset system in Eq. (58) the term x 

can be identified as the reaction torque due to the control, while 
Nq is the gravity torque. The magnitude of the reaction torque for 
the actual Z boom motion in the presence of gravity torque is compared 
with the magnitude of the gravity torque in Fig. 17. f. 

From these responses it can be concluded that the gravity torque 
has a profound influence on the system behavior. If the controller senses 
the actual angular velocities with the gravity torque present the very 
high amplitudes of the Z boom motion would force the spacecraft to a 
larger amplitude nutatlonal state. If the idealized control is used 
even though the response is a little better than the previous case the 
results are clearly undesirable. These results Indicate that with 
gravity-gradient present this method of control could be used for 
only a very limited time period, after which the Z boom motion would 
have to be terminated. 

It was thought that the presence of a counter mass, such that the 
composite center of mass would coincide with the hub center of mass for 
the zero position of the Z boom, might Improve the system response. 

To study this configuration the X boom mass, m 2 , was chosen equivalent 
to the Z boom mass and placed so as to maintain the symmetry of the 
spacecraft initially when z = 0 (see Fig. 16. a). 
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Identical initial conditions and satellite parameters as used in 
Fig. 17 were employed. The Z boom motion was assumed to be the 
ideal motion based on the control law for the torque free system. 

The Z boom motion and the nutation angle response are shown 
in Fig. 18. a. It is observed that with gravity torque even in the 
presence of a counter mass the satellite would achieve a very 
large nutation angle reaching as high as 80 degrees after 39 
seconds. When the gravity torque is neglected the nutation angle 
decreases initially and has almost a constant value of 2 degrees • 
after 20 seconds. Fig. 18. b shows the response of the components 
of angular velocity. The decrease in the magnitude of the 
component and increases in the magnitudes of the and com- 
ponents explain the high value of nutation angle at 39 seconds 
observed with the gravity torque. The responses without the 
gravity torque show a periodic motion for and after the 
initial transient, while remains nearly constant. This example 
with a counter mass points out the undesirable response in the 
presence of gravity gradient for the choice of parameters and 
initial conditions used. 

It has been pointed out in Ref. 7 that there is a residual 
oscillation in the Z boom after the spacecraft has been detuabled 
iising the control sequence. To study the response of the system 
in the presence of gravity torque, when the Z boom executes a sinu^ 
soidal motion the following equation is chosen to represent the 




I 


}, 






motion, of the Z boom: 


z = 4 sin 0.314t 

The responses with and without the gravity torqtie are shown in 
Fig. 19.a,b. In Fig, 19. a we see the Z boom motion and the 
corresponding nutation angle response. The initial conditions are 
indicated. The nutation angle response with the gravity torque 
shows periodic peaks reaching a maximum value of 9 degrees , . while 
intermittently reaching a minimum value of about 2 degrees. When 
the gravity torque is absent the nutation angle is periodic over a 
mpan value of 5.5 degrees. The corresponding responses of the 
angular velocity components are shown in Fig. 19. b. In the presence 
of the gravity torque and 0)2 intermittent peaks and troughs 

while Uo is oscillatory over a mean value of 0,55 rad/sec. The 

j 

responses without the gravity torque show that and are 
sinusoidal while remains a constant. The adverse effect of the 
gravity torque during the boom motion is clearly seen. 

A close study of Eq. (47) explains the substantial effect of 
the gravity torque observed in all the cases described above. The 
presence of the term results in the gravity force on each 

control mass exerting a torque about the hub mass center. Clearly 
the torque resulting from this is an order of magnitude in r/R 
larger, than that which would result if the torque were taken 
about the composite center of mass. Here it would be inappropriate 
to call the torqije as gravity-gradient torque since it is usually 



understood to mean the torque about the satellite composite mass 


center due to differential gravity force acting on each element 
of the body. 

In Eq. (47) if we isolate the term 


r 



X r. 


we see that even in the presence of the counter mass the vector 
addition of r^ and r^ (if r^ represents the position of Che Z boom 
mass and r^ the position of the counter mass) for the maximum z 
displacement would be such as to make the above mentioned term 
quite appreciable. This explains the undesirable response in the 
presence of the counter mass. 

From the previous cases studied it was concluded, that to Improve 
the performance of the system in the environment of the earths' gravity 
either the control mass and/or the amplitude of the Z boom oscillations 
have to be smaller. Further, the Initial tendency of the nutation angle 
to decrease even in the presence of the gravity torque suggests that 
the presence of suitable damping in the boom mechanism would serve to 
diminish the undesirable large amplitude residual oscillations of the 
Z boom, while maintaining Che initial beneficial effect of the control. 
Consequently two changes were made in the satellite parameters listed 
on page 36; the new mass of the satellite M = 4258 slugs (6.21x10^ kg) 
and the new Z boom control mass m = 55.95 slugs (816 kg) . Introduction 
of damping in the boom motion would result in the control Equation (64) 
being modified as: 




1 . 


2 + c (0.314) 2 + (0.314) 2 = -0.314 x 17.72 {kcj^Mj^+kc2t02J where 

c is the damping constant. The constants kc^^ and kc 2 can he determined 
for the modified satellite parameters as suggested in Ref. 7. They 
were calculated to be: 

kcj^ = -3.3188, kc2 “ 31.366. 

The response of the system for this case has been plotted in 
Figs. 20. a and 20. b under the same set of initial conditions and boom 
offset coordinates as used in Figs. 17. Figure 20. a shows Z boom motion 
for three different cases. The Z boom motion is indicated by the solid 
line in the absence of gravity torque and damping; by the dash-dot line 
in the presence of the gravity torque but no damping; by the dash-cross 
line in the presence of the gravity torque and in the presence of boom 
damping. The damping constant for this example was chosen to be, c = 0.4, 
which represents less than critical damping in the absence of control. 

The corresponding response of the nutation angle is shown in Fig, 20. b. 

In the absence of gravity torque and damping the nutation angle approaches 
2 ero within 175 sec. With damping in the boom motion the improvement in 
the response of the system in the presence of gravity torque is clearly 
seen. It should be mentioned here that the damping constant, c, must be 
chosen carefully to obtain the best results. 

Figures 21. a and 21. b show the response of the system for sinusoidal 



Z boom motion with the reduced mass. A comparison of this with the 
responses shown in Figs. 19. a and 19. b clearly points out the diminished 
effect of the gravity torque when smaller masses are employed. 


IV, CONCLUSIONS 


For deployment maneuvers during which the inertia symmetry 
about the spin axis is maintained the stability chart for 
symmetrical spinning bodies can be xised to study the 
system behavior. 

For various extension, maneuvers the bounds for boom 
lengths can be determined so that the system does not 
enter an unstable region. 

For extension along all three axes with the same end 
masses and extension rates, if the satellite is initially 
stable, there is negligible difference in the response of 
the system with gravity-gradient present or absent. 

For all the cases where the extension maneuver is performed 
through an unstable region although the nutation angle 
does not increase as soon as the system enters the unstable 
region, there is an exponential type growth in the nutation 
angle after a certain length of boom has been extended. 

The analytical solution obtained for the out of plane 
Euler angles for symmetric extension in torque free space 
can be used to study the system behavior with gravity 
torque for a limited time period. 

For the asymmetric deployment the expression for the 
gravity torque de\’eloped shows that a first order gravity 
force on the as3nmnetric mass exerts a torque about the 


hub mass center. This torque is an order of magnitude in 
(r/R) larger than the gravity- gradient torque for the 
S 3 mmietrlc deplojnaent. 

7. For Smaller control masses the controlled (asymmetric) offset 
boom system may be used effectively in the presence of the 
gravity- torque to reduce nutational motion when there is 
damping present in the offset boom system. The amount of boom 
damping required must be carefully selected for best results. 

8. For large control masses with gravity torque present the offset 
(asymmetric) system can be controlled so as to reduce the nutation 
angle only for the first few seconds, suggesting that the control 
should be terminated after that. 

9. The effect of other perturbations such as due to solar radiation 
pressure, aerod3mamic effects etc., have not been considered 
here but should be investigated especially for the case of large 
boom lengths. 
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DYC3) = C^2*i'il + ( 11-12) -DI3*w3+K8*CIl -12 )*CU2*C02*GD3*S053 /I3 
DY ta) = (i‘U*C03-W2*SD3 + Z*COV*3a2) /ca2 

. DY C5) = C'Yl+S03+if<2*Cq3-Z+S0U ... , 

DY t6) = (w3*C02-v<‘l*C03*S02+’fl2+SU2fcS03-Z*eui)/C02 
2Z=3QRrc ClI+ivn**2+(12*W2) **2) ■ 

^Y = I3*«3 .\,7 ■ 

ALFA = ATAi\i2CZZ,ZY)*57.2953 

HfcTUR'M 


PHaGRAMIS R£LOCATABLt . .. .. 

, TITL RGSOl . 
iFORT/A/H/E/P/5 FURl.LS/L 

TLISTING. . , .: , .. : 

7 SUBWtruriNE PGSrj2r.T,Y,DY,IHLF,N,P) 

7 LOGIC AL • 

; CUMMON .^n,Y-2.,y‘.3, A1,.A2, A3, ALFA 

; CUViMOM 110,120, 130, 1 1,12, 13,Ml, M2, M3rCl,C?,C3 

7 DlMt^JSla^i YC6),DYC6),0UMMYCb) 

.7...C.,.:. .... . ■.. . . : V- ; . ■ ^ ■ .■ . . .- 

7 CALL RGSOl CT,Y, DUMMY) 

7 IH=T 

7 .. . .lFC*'^OT.P.Kr>.XUlHLF J) GU TO S . . 

: lYK 1 T £ C S , 1) TP , Y c n , Y C2 1 , Y ( 3 ) , Y C a ) , Y c 5 ) , Y c b ) , AL F A , IHLF 

7 1 FURMATClx,F7.3,btlX,El3,8),lX,Fl0.7rlX,I3) PAGl 3^ 

7 . HlMAPY (1) T,YCn,YC2),ALFA . OF POOE QUAIM 

7 8 CO'ir IMUt 


I OOOOOOO OO * 1 3ZIS 
; 000000000*0 A 

OOOOOOOSjO*0 000000090*0 000015000*0 000090000*0 000090000 * 0 9ZTS 

U00000000*0 000000100*0 000905000*0 000000000*0 000050000*0 A 

000000*0 =£3 000000*0 =23 000000*0 =13 

000000*01 =0£I 000000*5 =021 000000*9 =OlI 

0000T0*0 =£w 000000*0 000000* 0 =TW 

• ■ , ‘ 3'9X'^T 

8TJH09 ^"I'cJSStOda £00 200 lOD S/dlMi I'O/HcnHt 
2[JS3a liii* 

939VJLV3n‘13H ST WuasnHd 


■■ '• ' ■■■ ' ; ■ " ■■■■"■ . ' 
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ASmiETRIC DEPLOYMENT 


I5iba;22 

i JOB, . tRg AD..J N At.: .if ? ; 53. j S.9J 
IFORT/6/b/e/P/S FOPT.LS/L 
iL.TSTIf.‘G 




A5Y‘'*»‘'t:TffIC. APPENOAGES 
RG=1 .0 


C . EFFECT np CPAVHY TORQUE 
C ALT1U.0F .= 500 iiAUT.lCAL MILES. 

C INCLUnT.VG GRAVITY-C-RAOIEMT EFFECTS 

EXTEP-'J.AL . R.GSUl ,PGSri2„. '. .. _ 

OXME^i■SIO^^ PARV(51 ,y (6) ,DY(fa),waPKCa,6D^SI7E (b!) 

REAL 

. . CO‘-)^OA! Y . ■ 

CUM^^o^J LI, 12,13/ '-'S> Ml , M2 

COMMON A A, BS, CCrDO 

com miim. Z.07. nn7 : 

COMMON N, M 
COMMON CC3/3D 

..... .. ..... .E<3UIv.ALEiNCE.CY..C.U/.w 1.1/ 1 Y C 2) / W 2 ) ,. CY.C3 ) , W3 ) /.{ Y (« ) , A n / CY.C 5.J , A2 J 

3CYC6) , A3) 

CALL INCUJTC2/5) 

C.ALL...QP£N.a.,„': R A JAN .3/ 1 ER ) . 

IFCIER.NE.DSrOP LjrviflRLE TO OPEN FILE 
HEA0(2/9UTMA.X,STEP/T0L. 

9i._ ..format c.aFi 0 .0) ;. ....... ; .. ... : .. 

PARM(i)r0,0 ■ 

PARM(2)=TMAy : 

pAjiNi,aiF.siEP_ 

C INITIAL VALUES 

R£40C2,91)Y . • 

.. :.. .--....read C 2/.R 1 ) SIZE.:. . : ■ 1 ; ... : .. . 1 ... : ............. ........ . 

PEA0C2/PU)I1/I2/I3 

REAOC2/9t nM8,Nn/M2 ‘ 

...... RE.AD.C 2 ,..9 lilA.A , SB /.CQ.,.DO .......... 

911 FORMATCOFPO.A) 

N=b' ■ '1. ^ ‘ ' 

- 1 .... . 1 . 

wRiTE(5,9 21TN'Ay /STEP, TOL 
i'siRlTEC5,93)l3 ,12, 13 

■_ 

i'iRITFcS/95)AA,8B,CC,D0 
" rtRlTECS,9b)Y 

'..WRITE C5,97)SII£........ ... , I' . ....I.- 

ivKITE (5, 98) 

9? FOR'-'' AT CIX , 'TYAxs' , F 1 0 . b , 5X , * ST EP= ' , FI 0 . 6 , 5X , 'TOL = * ,F1,0 ,6) 

. : ..9 3.. .._.FOi?3 A m X , .'.U F /.E 15.*. 4. , 5 X ,..U 2= ' ,.£15 . Xj,.: 13= '. /.F-.l 5 . 0 ),_. 

9U FOR'VATCIX, 'mR = ',E1S. a,5X, 'f-n = ' /E15.9/SX, *M2E:',El5.ii) . 

. 95 F OR M A Y f I X , ' A A= ' , F I 0 . 6 , 2X , ' 8B= ' , F 1 0 , 6 , 2X , ' C C= ' , F 1 0 . 6, 2X , ' DD= ' 

i£'Fl0...6l - . , . .... 

9b FGRMATC3X, 'Y ',3X,6F15.6) 

97 FURMATC3X, 'S.IZE' ,bF15.6) 

.. . .9 S ..:.F.Q.9ZlA_t C .' 1.' , TB, .1.11, T.L.7 ,_' ill .1, 13,0 / l.Y 2 1/.TA3 , ' iv3.' ,.T 55 JHET A.- 1.^. 
2T69, '1 HETA-'2' ,T8 1 , 'THETA-3' , T9 6, ' ALFA ', Til [) , ‘’IHLF' ,/) 

C 

CIALL.._9»vSCL.tiV , S.I 7F., UY , T.OL, P-ARM) ;... , 

CALL R.MGSCPARf--., Y , DY,N, IriLF,RGS01 , 90502, wpPX) ' . i, - ■.'.'I 

i‘‘9ITEC5,99)IHLF 

.....99.. F09,MATC'f)I.HlF= 1,1 3 )....,. .l.-.l 

Call exit 


'.-ij 


PfiOG(^A^‘ IS wEL0CATA6Ue 

,T1TL , 

iFORT/ft/K/e/P/S F0PT.L3/L • 

ILISTIHG 

SUARnunwf RGS_01 IT, Y^OY> . . _ . 

L>ivE(^^SIO^: Y Cfr) #DY ChD rCC]^,3) 

PEAL 

. COM.viuw :/1,A2,W3/At/A2,A3 ,• . . 

CU-V^viON IV,I2rl3,Me,Ml,M-2 
CUM:^)n^} AA,«R,CC,00 

Cr^vlViONj ZrPZ,,P{>2 . ,... 

CO.'^'''iOM >I,P' ■ 

C 

c... . CAL. ..OF CQEFF..-.OF L.ri,S. . OE MATRIX EpM, 

Ul=Ml*C--^6 + M2)/{MB+Mt+M2) 
g2=M2*CMB+Mn / (W8+M1+M2) 

.U3.= -M_1*M2/C.MR.+Mi+M2). .. . _ . _ _ _ 

X=O , 0. ' . 

DX=d.Q 

DDXgQ^n „ 

TTi=0. 314*7 

S(sM=SIN(TT) 

C.O A =C..Q.SXT.T 7 _ _ ■ . _ : ... _. 

er = EXP(-0.5il*TT) 

CO2=C0S(0,a28*TT) 

ii.-L.^:._sy2.=s.iiiC.o»^ ..... ... 

Qi = 0.655*s'Ql + 0.a{7*CQl 
O2=ET*(0.al7*CQ2+?.056*SQ2) 

.Q3.=£X*CrQ..AIl*Q..A28*SG2+2.J.56*0,4.2ft*CQ2). _ „ _ 

132 = 01.-02 
2=17.72^^02 ■ . 

oDZ,?n-.(>.55*coi-o..4i.7.*so.ir03+o...5A*Q.2 _ _ 

02=0.314*17, 72*QDZ 

QDOZ=-Ql-ET*C-0.iil7*0,A2S*0,428*CQ2-2. 056* 0.428*0,428*8627+2.0* 

®7X5.4_*9.3.-.0., gt|*0.Ba*Q2 

DD2=Q.314*0.314*17.72*QODZ 

CCl, n=IlfUl*C8B**2+2**23+U2*(CC**2+DD**23t2.*U3*C8B*CC+0D*Z) 

- . C Cl . , 27.=-U 1*A A.*BS-U2*.CC.*X.-U 3* ( A A*CC;+8fi*X) . . _ . ■ 

Cllr31=-Ul*AA*Z-U2*DD+X-U5*CAA*DD+X*Z) 

C(2, n=GCl,2) 

-__L.___C.C2,X7_=I2tjUL*CAA**2+Z**2 7+t)2* CDD**2*X**21*2.0*U3*CAft*X+DO*Z7 

CC2f 3)=“Ul*8B*Z-U2*CC*D0-U3*CRH*DntCC*Z) 

CC3> n=CCi,3) 

CC3./21TCC 2, 33.... . ........ . . , .... 

C i;3»3} = I3+Ul * C AA**2+BB**2) *U2* CCC-**2*X**23 +2,0*U3*C8B*CC + AA*X ) 

C CAL. OF P.H.S OF MATRIX EQUATION 

..L..... .. .:GT=l,A07528:*.tm*Ji**^^^ .. . .. . 

R=2.3P508*C1 0,G**7 ) ' 

: GG=Un. . . 

. ^=GT*GCi : : ... .. :.... . 

ii.'«=5QkT(GT/(P**3) ) 

cai=cns(Au _ 

_.„__C0 2=008 C.A 23; 

c05=CnsCA33 PCK 3 B QlJAIiM 

SOI =STfiCAl 3 • ' 

. 8U2=6.INCA2) ..... ... ,. . . ; 

SU3=STnC A3) 

P*< = 3,0*K/P+*3 

H1K=Pk*m1 * C.a a.*C02*C.03.-R8*002*$03*2*5023 , ... . .1 ; . 

R2><:=p!l*'-^2*tX*Cn2*C0 3-CC* 002*505*00*5023 
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f . 

^ . . 

w 

f 

* 

r 

« 

r 


Z0m=2*mi +OD*-Ma 
KC«st^R*M1 +cc*^^a * 

. Ai l=UJ ^C-AAitZ-*'/*! *VV2+ teR*'t2-Z**2.)*w2*V'3+AA*BB*Wi*W3 + 2,*2*0Z*Wlt .. 
*bB*Z*fv^!5*-*2-''/2**? ) +Hrt*OOZ) 

A12=:U2* C-l.'>D*X*i‘ 1 *?^2 + CCC*>2-D0**2) *V'<2*;h3+CC*X*Wt*W3-2,*CC*DX*i*i2 

. .. . ■«-2.A06*0X*«?.3+q C*. D P *.. ( W 3 * * ? - W ? *-iit 2 1) 

A1 5=U3*( C-AA*DD-Z*X3*vd;ifVV2 + 2,*((5B*eC-D0+Z)*W2*iM3+CAA*CC + Be*X) 
*«3+?,*nD*DZ*rti-2,*8B*L'»X*’/J2-2,*DX*Z*irV3+(8B*00 + CC*Z)*Cv'J3**2 

+CC*(>OZ) . 

Al £|=HK*( l?-l33*Ca?*SQ2*S03 
A1B=CK/H**23*C-C02*S0 3*ZD’M-S02*6CM) 

. . A16=RlK*CZ*C0E+S03+-Ba*SQ23 . , 

A17=R2K+(OD*C02*S03+GC*SQ2) 

l)y t 1)-- CA1 OtAi 1+A12+A13-AU-A15-A16-A173 

BlD-CIl-I3)*;-v3*V'-l, . . 

C-BB^ZtHnI *w2 + AA*Bt?*v^2Vw3+( AA**2-Z**2)*W3*vn-2.*Z*DZ*W2 
-i)+AA*Z*f w3#-*2-Kl**2) +AA*OOZ) 

.8 1B=1J2.*Cj- CC..>r.Q£i t.l'v.l* l^g ■fcC.C.*.X »W2.*W3.+_C X **2-D0^* 2 ) *W 3*tN 1 -g . * X*DX * W3 

W+DX*(i''J5**-2-*'<i**2)”OD*DDX ) . 

013=U3* (-CBB*DD+CC*Z 3 *v^ i */J2+ ( A A*CC+BB*X3 *W2*W3t2 .* C A A* X«DD*Z) 

3*vy.3 *>'( L-Z..M A A. * 0X>D.P *aZl* « 2.t,.C A A *:QD.+ X *Z )jLCJ« 3.*^ 2r W 1 **Z') + X..* D D Z - Z * D 0 X ) 
8ia=RK* Cri-13D*S02*C02*C03 

bl5=CK/R**23*(-C02*C05*Z0M+S02^cAXM) 

B ib~ H iKAc z* CO 2 .*c 03 r A A*.aa 2 j : 

B17=R2K*CDD*C02*CQ3“X*C02) 

DYC2)--CB10-eil»ai2-at3-B14-815'^816-B17) 

Cl 0 = (12-1 U*)'- 1 m2... ..... 

Cl l = Ul*' ( C86**2-AA**2}*Wt*w2’-AA*Z*W2*W3 + B8*Z*W3*Wl + 2,*AA*D'z*WH> 

dt2,*B8*-DZ*!,N2+AA*88* CWl7(f*2-w2**^23 ) 

... Cjt ?=>J2*C fCC^*2-X**23 *l^l*'^2-DD *X*tM 2*vM3-t-CC*DD*W3*Wl-2 .*X*0X*rt3 

n) + C C * X ■*■ ( 4'-' 1 * * 2 — 2 * ■* 2 3 + C C * D 0 X ) 

C13=U3* C2.*C8B*CC-AA*X3 *wl*W2-(AA*DD + X*Z3*W2*W3+(a8*OD+CC*Z3*W3*'Wl 
^1-2. A CX*nZ*Wl +-CC*DZ*W2-AA,>0.X*W3).±..(AA*CC+B.fii»rX)*.(Wl*^.2-W^^ -f-aB* „„ 
SDDX3 

CiqsRK* Cl 1-123 *CD2*C02*C03*S03 


C 15.= C C.0a*SJl3..*MM*Cl)2ia:Q3.*.ec.Mi; 

j C16=Rj K*CAA*C02*Sa3+BB*C02*C033 

; ■ CI7=R2K*CCC*C0P*GO3+X*eO2*SO33 

? 0YC33.=-CC10-CUrCl2-C.1.3«C.14-CL5'trC.lb + Cl7) .. .. 

; UY(4) = C'?Jt*C03-t'J2*S03+0M*-C01*S02)/C02 

; DYC5) = Cv-il*S03+i^j2*C03-0M*S0n 

; ..D Y.Ld,).= C4'J3*Ga2r3'^i * C a3*.5 0.2 4-A 2*SU2*503-0M*CQ 1 3 /C02 

; CALL SIMOCCrOY,N,KS) 

; IFCKS33,273 

2 .. RE1.U.RM ... 

? 3 wRITECBf 4) 

; £1 FO«fflATC//*SIKGULAR EQUATIOMS'3 

; . . . _ . R EI.U B N_.: . .. , : , ._ .. . . . .. .:.... 

; - £N0. 


PROGRAM IS relocatable 

.TITL RGSOl 
*P0 RT/a/p/£/P/. 5. PQRT,L5/L 
iLISTlWG 

: SUPBOUTIRE FGSO?CTrY/OY, IHLF,N,P3 

: . ... . . LOGICAL f-KGX T . . . . ...... 

; 0 1 ^lE MSI OR Y C A 3 , F3 Y C b 3 , DUMMY ( 61 

; HP^AL I U12,I3/f-BfMlr-’2,K 

5 . COMMON . M L I w2# 03 r A I t A 2. r.A3 

: II Tl2Tl5f <^’HT ^'1/1^2 


PAG® ® 


0g ^OOB 



CflM^OM fi/i , KK,CC , DO 


; C 

j CALL kGSIJ] CT, Y jOUV’-iy ) 

j tHl = Il*-lNj, 

? H?=I 2 *'A ;2 

; H 3 = 1 ? * '■i ? 

; A LF A = fiWPf SORT CHI ,H3) *57,2958 

; TP = r 

; ^FC,l■•KJ^.RK^;!<TCTHLr n GO TO 8 

; 1 ) rP,Y ,7,07, ALFA, IHLF 

; 1 KS>Wvi/^T(;iX,F5,l,8nx,F12.5),lX,F7.«,lX,I?.J 

; WRITE 51MARY (1) T,Z 

; 8 CUV r I WUE 

; RETURN 

: END 


PROGRAM IS REL0CATA6LE 

,TITL RGS02 

iRLOR/M rwP/S Ooi 002 003 PPOjSSP.LB F0RT.L8 
‘EXEC 

TMAX= 60.000000 STEP* 3,000.000 T.OL= 0.010000 

li=. 0.1050E B 12= 0.1050E 8 13= 0.1500E 8 

MB= . 0.1570E 6 Ml= 0,18006 0 M2= O.OQOOE 0 

AA.= .65,J)On_0.00 ..BP= 0,000000 CC= 0,0000.00, OD= 0.000000 

Y 0. 039100 0,000000 0,31^000 0,100000 

SIZE 1. 500000 1.50000Q 2,000000 2,200001 


Y o'.oooobb oV'oooooo 

SIZE 1.000000 .42.000000 



obkhnai) f ® 
6 E EOOE (JUAIiT^ 




